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Abstract 


This paper deals with the estimation of P(Y < X), where Y has generalized 
exponential distribution with parameters a and A and X has mixture gen- 
eralized exponential distribution (or marginal distribution of X1, X2,...,Xn, 
in presence of one outlier with parameters 6; and $2) such that X and Y 
are independent. when the scale parameter (A) is known the maximum like- 
lihood estimator of R= P(Y < X) is derived. Analysis of a simulated data 
set has also been presented for illustrative purposes. 
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1 Introduction 


Recently the two-parameter generalized exponential (GE) distribution has been 
proposed by many authors. It has been studied extensively by Gupta and Kundu 
({11}-{17]), Raqab [26], Raqab and Ahsanullah [27], Zheng [34] and Kundu et 
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al. [21]. Note that the generalized exponential distribution is a submodel of the 
exponentiated weibull distribution introduced by Mudholkar and Srivastava [22] 
and later studied by Mudholkar et al. [24] and Mudholkar and Huston [23]. 


The two-parameter GE distribution has the following density function 
fede oe “lee "4, ie Oe (1) 
and the distribution function 
F(az,a,A) = (1 —e7**)?, z>0 (2) 


Here a > 0 and » > 0 are the shape and scale parameters, respectively. For 
different values of the shape parameter, the density function can take different 
shape. For detail description of the distribution, one is referred to the original 
paper of Gupta and Kundu [11]. From now one GE distribution with the shape 
parameter a and scale parameter 2 will be denoted by GE(a, A). 

Let the random variables X 1, X2,...,Xn—1 are independent, each having the 


probability density function f(x), 


filz, Bz) = Poe *(1—e *)P1, o> 0, (3) 


and the one random variable (As outlier) is also independent, has the probability 


density function g(x). 


folz, B1) = Bye *(l—e *)t, z>0 (4) 


The joint density of X1, X2,...,Xn is given as 


= Tp Fey. y> Bead 
7 (Ei toe) = au ea VU ia) 


— 

i ie ae Bye *41 (1 — e241) 1-1 
= er — v1 Xi = e *) B2-1 1 

me e IIa > ( Boe” TA, (1 —e *41)62 =] 

1 rm nm 
_ = By Bo" 1e7 wey i [a = e ti)Pal > (2 e741 )P1—Ba 

n 

i=1 Aj=1 


,c>0 (5) 
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(see [7] and [8]) 


From equation (5), the marginal distribution of X is, 


n—1 


eae <, «>0 (6) 


h(x, Br, Ba) = Bre") 


n 


The main purpose of this paper is to focus on the inference of R = P(Y < 
X), where Y ~ GE(a,A) and X has mixture GE or marginal distribution of 
X 1, X2,..., Xn with presence of one outlier. For simplify we consider 4 = 1. The 
estimation of R is very common in the statistical literature. For example, if X is 
the strength of a component which is subject to a stress Y, when R is a measure 
of system performance and arises in the context of mechanical reliability of a 
system. We obtain the maximum likelihood estimator (MLE) of R. It may be 
mentioned here that related problems have been widely used in the statistical 
literature. The MLE of P(Y < X), when X and Y have bivariate exponential 
distribution, has been considered by Awad et al. [2]. Church and Harris [4], 
Downtown [6], Govidarajulu [9], Woodward and Kelley [33] and Owen, Craswell 
and Hanson [25] considered the estimation of P(Y < X), when X and Y are 
normally distributed. Similar problem for the multivariate normal distribution 
has been considered by Gupta and Gupta [10]. Kelley, Kelley and Schucany [18], 
Sathe and Shah [29], Tong [31], [32] considered the estimation of P(Y < X) when 
X and Y are independent exponential random variables. Constantine and Karson 
[5] considered the estimation of P(Y < X), when X and Y are independent 
gamma random variables. Sathe and Dixit [8] have been estimate of P(Y < X) 
in the negative binomial distribution. Ahmad e¢ al. [1] and Surles and Padgett 
[30] considered the estimation of P(Y < X), where X and Y are Burr Type 
random variables. Baklizi and Dayyeh [3] have done shrinkage estimation of 
P(Y < X) in exponential case. 

The rest of the paper is organized as follows. In section 2, we derive the MLE 
of R. Analysis of a real life data set has been presented in section 3 and finally 


we draw conclusion in section 4. 
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2 Maximum likelihood estimator of R 


Let Yj, Yo,..., Ym be a random sample for Y with pdf 
swe) Soe Kise 4s y>0 (7) 
and X,, X9,..., Xn, be random sample for X with pdf 


n—1 


Bet See, 0 8) 


f(x, B1, B2) — je*(1 _ Cae ae a 


Then 


n 


he => PW) 


n 


+ i [ec — | pe-*(1 — 6 *) "ae 4 (1 — e749" ]dydax 


nat By n at Bo 
Therefore, the MLE of R becomes 


: 1 6 -1 £6 
Rk = eae eS ae 
Nat By nm a@+t Bo 
Now, to compute the MLE of R, we first consider the joint distribution of 


X 1, X9,..., X,with presence of one outlier in (5), so 


L(a, Bi, B2) GS Uo eg Ue) Ege OG) 
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The Log-likelihood function of the observed sample is 
m m 

In L(a, 81,82) = mln(a) — Sout (a — 1) S$ n( —e %) 
i=1 i=1 


ipa * n ~ 
+ In[——e- wisi me (= e ti) P2-1 S> ie 7 ti) P1— 82) 
Ai=1 
(10) 
The MLE’s of a, 6, and (2 say a, By and Bo, respectively, which is obtained as 


the solutions of 


OlnZL m = 
— ] 1-— Yi) — 
- 2 nile") 0 


Oa 
or 
m m 
ras ee —e %) 
ra i=l 
Hence 
—m 
A a 11 
“= ST In(l — e-¥) (11) 
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Form (13), (14), 8, and fy» can be obtained as the solution of non-linear equations. 
For 2, = 2 = P in homogenous case, @ and B can be obtained as 


as —m B _ —n 
© SP In(l =e) ~ Se In(l —e-¥) 


These equations proposed by Kundu and Gupta (2006). 


(14) 
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3 Numerical experiments and discussions 


In order to have some idea about Bias and Mean Square Error (MSE) of MLE, 
we perform sampling experiments using Maple. 


To have inference about R, we consider the following small sample sizes; (m,n)= 
(15, 15), (20, 20), (25, 25), (15, 20), (20, 15), (15, 25), (25, 15), (20, 25), (25, 20). 


Here, we take a = 1.50 and 6, = 2.5 and $2 = 2.75, respectively. As we know, 
the generated sample size n from f(x, 61,2), (n-1) sample generated from the 
equation (3) and one sample generated from the equation (4). All the results are 
based on 1000 replications. Here we present a complete analysis of a simulated 
data. The data has been generated using m = n = 20, a = 1.5, 6, = 2.5 and 
Bo = 2.75. 

The data has been truncated after two decimal places and it has been pre- 


sented below. The Y values are, 


0.74 1.41 0.86 0.20 0.72 3.11 0.73 0.44 
1.31 0.86 0.27 2.27 0.88 1.32 4.41 1.17 
0.86 2.19 0.53 0.08 


and the corresponding X values are, 


1.12 5.30 0.65 1.46 1.27 0.74 1.51 0.81 
1.79 2.11 1.33 1.50 1.57 1.26 0.49 2.93 
0.85 0.85 1.73 1.83 


Now, we obtain the MLE of & = 1.671, 6; = 0.315 and 8) = 2.18. Therefore, 
R = 0.5457. 
4 Conclusions 


In this paper, we have addressed the problem of estimating P(Y < X) for the 


Generalized Exponential distribution with presence of one outlier, when the scale 
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parameter is know. The results are given in table 1, 2 and 3. It is observed 
that the maximum likelihood estimator R work quit well. We report the average 
estimates and the MSEs based on 1000 replications. The results are reported in 
the following Tables. In this case, as expected when m=n and m, n increase then 
the average biases and the MSEs decrease . For fixed m as n increases the MSEs 


decrease and also for fixed n as m increases the MESs decreases. 


Table 1 
a= 1.5, By _ 2.5, Bo = 2.75 


a 

(15,15) 0.1837793490 
(20,20) 0.1882900337 
(25,25) 0.1811570295 
(15,20) 0.1767666919 


(20,15) 0.1804719689 
(15,25) 0.1719679412 
(25,15) 0.1791734199 
(20,25) 0.1854977684 
(25,20) 0.1847561247 
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Table 2 
a= 1.5, By = 2.5, Bo = 2.75 


oa or 
(15,15) 0.2526986879 
(20,20) 0.1554931754 
(25,25) 0.1135263151 
(15,20) 0.1467095745 


(20,15) 0.2182551913 
(15,25) 0.1069492224 
(25,15) 0.2132163252 
(20,25) 0.1176980818 
(25,20) 0.1455059822 


Table 3 
a= 1.5, By _ 2.5, Bo = 2.75 


oi = 

(15,15) 383.2417437 
(20,20) 14.25371765 
(25,25) 1509.645815 
(15,20) 19357.50094 


(20,15) 20.41377913 
(15,25) 33.74240776 
(25,15) 628.5131946 
(20,25) 86704.72425 
(25,20) 5162.227553 
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Table 4 
a= 1.5, By = 2.5, Bo = 2.75 


(a.m) 
i615 
(20,20) 0.6090520559 
(25,25) 0.4601051183 
(15,20) 0.9969076983 


(20,15) 0.5981276207 
(15,25) 0.8879485964 
(25,15) 0.4227174191 
(20,25) 0.6328966532 
(25,20) 0.4015675626 
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